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Abstract: Eco-environmental quality is a measure of the suitability of the ecological environment for 
human survival and socioeconomic development. Understanding the spatial-temporal distribution and 
variation trend of eco-environmental quality is essential for environmental protection and ecological 
balance. The remote sensing ecological index (RSET) can quickly and objectively quantify 
eco-envitonmental quality and has been extensively utilized in regional ecological environment assessment. 
In this paper, Moderate Resolution Imaging Spectroradiometer (MODIS) images during the growing 
period (July-September) from 2000 to 2020 were obtained from the Google Earth Engine (GEE) 
platform to calculate the RSEI in the three northern regions of China (the Three-North region). The 
Theil-Sen median trend method combined with the Mann-Kendall test was used to analyze the 
spatial-temporal variation trend of eco-environmental quality, and the Hurst exponent and the Theil-Sen 
median trend were superimposed to predict the future evolution trend of eco-environmental quality. In 
addition, ten variables from two categories of natural and anthropogenic factors were analyzed to 
determine the drivers of the spatial differentiation of eco-envitonmental quality by the geographical 
detector. The results showed that from 2000 to 2020, the RSEI in the Three-North region exhibited 
obvious regional characteristics: the RSEI values in Northwest China were generally between 0.2 and 0.4; 
the RSEI values in North China gradually increased from north to south, ranging from 0.2 to 0.8; and the 
RSEI values in Northeast China were mostly above 0.6. The average RSEI value in the Three-North 
region increased at an average growth rate of 0.0016/a, showing the spatial distribution characteristics of 
overall improvement and local degradation in eco-environmental quality, of which the areas with 
improved, basically stable and degraded eco-environmental quality accounted for 65.39%, 26.82% and 
7.79% of the total study area, respectively. The Hurst exponent of the RSEI ranged from 0.20 to 0.76 and 
the future trend of eco-environmental quality was generally consistent with the trend over the past 21 
years. However, the areas exhibiting an improvement trend in eco-environmental quality mainly had weak 
persistence, and there was a possibility of degradation in eco-environmental quality without strengthening 
ecological protection. Average relative humidity, accumulated precipitation and land use type were the 
dominant factors driving the spatial distribution of eco-environmental quality in the Three-North region, 
and two-factor interaction also had a greater influence on eco-envitonmental quality than single factors. 
The explanatory power of meteorological factors on the spatial distribution of eco-environmental quality 
was stronger than that of topographic factors. The effect of anthropogenic factors (such as population 
density and land use type) on eco-environmental quality gradually increased over time. This study can 
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serve as a reference to protect the ecological environment in arid and semi-arid regions. 
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1 Introduction 


As climate change and human activities intensify, ecological environment problems have gradually 
threatened regional ecological security and socioeconomic sustainable development (Liao et al., 
2020). Together, the three northern regions of China (i.e., the Three-North region), including 
Northwest China, North China and Northeast China, constitute one of the most crucial natural 
resource reserves in China and form an important ecological barrier. With nearly 280,000 km? of 
suitable forestland (accounting for nearly two-thirds of the total suitable forestland in the 
country), the Three-North region has great potential to increase forest resources and is the key 
area for forestland and grassland construction and carbon sinks in China (Zhang et al., 2021a). 
However, the natural environment in the Three-North region is fragile and variable, and most of 
the region is located in arid and semi-arid climate zones (Zhang et al., 2021a). In 1978, the 
Three-North Shelterbelt Program, which covers 13 provinces, autonomous regions and 
municipalities in China, was launched to alleviate soil erosion and wind-blown sand hazards in 
the Three-North region (Duan et al., 2011; Deng et al., 2019). Since the implementation of the 
program, the ecological status and economic benefits of the project area have been greatly 
improved, playing an important role in maintaining national ecological security and promoting 
socioeconomic sustainable development (Cao et al., 2020). In recent years, to enhance the carbon 
sink potential of forestland and grassland in the Three-North region and achieve carbon peaking 
and carbon neutrality, the Chinese government has intensified efforts to protect and restore 
important ecosystems through the Three-North Shelterbelt Program (Lu et al., 2018); therefore, 
eco-environmental quality in the project area has been continuously improved. Eco-environmental 
quality is the comprehensive performance of ecosystem elements, structure and functional 
characteristics. Superior eco-environmental quality fundamentally guarantees human survival and 
is the physical prerequisite for socioeconomic progress (Jiang et al., 2021). Mastering the 
spatial-temporal distribution and variation trend of eco-environmental quality is helpful for 
ecological protection, environmental restoration and policy formulation in the Three-North 
region. It is important to comprehensively promote the carbon sink capacity of the ecosystem in 
the project area and to achieve both the 'double carbon' goal and sustainable development (Lu et 
al., 2018). 

The traditional evaluation of eco-environmental quality has primarily been based on the 
statistical yearbooks (Shan et al., 2019; Wang et al, 2019). Annual data were recorded 
systematically and continuously; however, basic data points regarding vegetation, soil and 
hydrology over large areas on monthly or quarterly scales have not been collected. Annual data 
are limited in their usefulness for the dynamic monitoring of eco-environmental quality over short 
intervals (Zhong et al., 2020). Remote sensing technology has the advantages of fast, real-time 
and widespread coverage. It also offers technical support for the monitoring and evaluation of 
eco-environmental quality (Zheng et al., 2020; Zhou and Liu, 2022). Monitoring indices 
constructed using remote sensing technology have become the primary means of the monitoring 
and evaluation of regional eco-environmental quality. For example, the normalized difference 
vegetation index (NDVI) has been widely used to evaluate vegetation growth and phenology 
(Erasmi et al., 2021; Zhang et al., 2021b), the water index has been applied to evaluate drought 
status or habitat suitability (Mishra and Pant, 2020; Al-Quraishi et al., 2021), and the land surface 
temperature (LST) has been used to evaluate the urban heat island effect (Lee and Park, 2020; 
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Arshad et al., 2022). The ecological environment is composed of complex ecosystems and is 
affected by multiple factors; therefore, a remote sensing evaluation method based on a single 
monitoring index cannot comprehensively capture the systematic changes in eco-environmental 
quality (Wen et al., 2019; Zhou and Liu, 2022). The Ministry of Ecology and Environment of the 
People's Republic of China proposed the ecological index (EI) to evaluate the regional 
eco-environmental quality in China from the five aspects of biology, vegetation, hydrology, land 
and pollution (Shan et al., 2019; Zhou and Liu, 2022). However, the EI can only provide a general 
overview of the ecological situation of a region, rather than a detailed description of the local 
distribution of various environmental conditions in the region. Based on the EI and remote 
sensing data, researchers constructed the remote sensing ecological index (RSEIJ) using principal 
component transformation to couple measures of greenness, humidity, dryness and heat (Hu and 
Xu, 2018; Xu et al., 2019). The RSEI can overcome the spatial visualization difficulties of the EI 
and reduce the impact of anthropogenic factors in the evaluation process (Hu and Xu, 2018; Xu et 
al., 2018). It has been widely used for monitoring eco-environmental quality in cities (Wen et al., 
2019; Sajjad et al., 2021), lakes (Xiong et al., 2021; Yuan et al., 2021) and other regions. 

The dynamic analysis of time series remote sensing images is important for identifying the 
changes and internal laws of eco-environmental quality. The existing research methods for 
long-term sequences are numerous and have been used frequently, providing reference for their 
application to the spatial-temporal analysis of eco-environmental quality (Sun et al., 2020; Li et 
al., 2021a; Li et al., 2021b; Xu et al., 2021; Jiang et al., 2022). The linear trend analysis method is 
commonly used to describe the trend and change of time series data by obtaining the slope to 
estimate the magnitude of the variation trend through regression analysis (Duan et al., 2011; Deng 
et al., 2019; Xu et al., 2021). However, this method has limited stability and is sensitive to 
anomalous data. The Theil-Sen median is a nonparametric statistical trend estimator that is often 
combined with the Mann-Kendall (MK) test for the trend analysis of long-term sequence data (Li 
et al., 2021a; Zhang et al., 2021b; Jiang et al., 2022). Compared with the linear trend analysis 
method, the Theil-Sen median trend method combined with the MK test has high calculation 
efficiency and robustness to error, which effectively reduces the influence of data outliers and 
improves the accuracy of test results. In addition, the Hurst exponent can be calculated to predict 
the future trend of time series data (Li et al., 2021a; Li et al., 2021b; Rivas-Tabares et al., 2021), 
and the coefficient of variation (Li et al., 2021b; Zhang et al., 2021b), Moran index (Xiong et al., 
2021) and change vector analysis (Xu et al., 2019; Sun et al., 2020) can be effectively used to 
understand the spatial evolutionary characteristics of vegetation. 

Natural factors control the evolution of the ecosystem environment, while human activities can 
have a considerable impact on the ecological environment of a local area. Therefore, a 
combination of natural factors (temperature, precipitation, etc.) and anthropogenic factors 
(population, land use, etc.) may contribute to the changes in eco-environmental quality (Gao et 
al., 2020; Sun et al., 2020). Correlation and partial correlation analyses are commonly used to 
determine the relationship between drivers and response variables (Hu et al., 2021; Liu et al., 
2021). Correlation analysis can assess the influence or association between two factors without 
considering the influence of other factors (Hu et al., 2021). Partial correlation analysis can 
identify the interaction between specific variables after controlling the influence of other factors 
(Liu et al., 2021). Correlation and partial correlation analyses are limited in their ability to 
quantify the coupling effects of multiple factors; however, the geographical detector can 
overcome this problem. The geographical detector, which is based on the factor force 
measurement index, can detect spatial differentiation by combining spatial superposition 
techniques and set theory to reveal the driving forces behind the differentiation mechanism (Wang 
et al., 2010; Wang and Hu, 2012). 

Optical data, such as Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat 
images, are commonly used data sources for constructing the RSEI (Shan et al., 2019; Wen et al., 
2019; Liao et al., 2020; Zheng et al., 2020). Due to the influence of weather, topography and other 
conditions, it is difficult to obtain cloudless Landsat images of a large area in the same period (Ji 
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et al., 2020). MODIS data have high spatial resolution and represent a complete time series, but 
they require excessive data preprocessing and index calculation processes when applied on large 
regional scales (Xu et al., 2019; Zheng et al., 2020). Google Earth Engine (GEE), a remote 
sensing cloud platform, contains massive numbers of historical images and geographic databases 
and has overcome the problems of low efficiencies of local downloading, storage and 
preprocessing (Gorelick et al., 2017). GEE has been applied extensively for large-scale regional 
data fusion (Nietupski et al., 2021), crop yield estimation (Adrian et al., 2021; Venkatappa et al., 
2021) and land cover dynamic analysis (Floreano and de Moraes, 2021). 

In this study, GEE was used to obtain and preprocess MODIS images during the growing 
period (July-September) from 2000 to 2020 to calculate the RSEI for the dynamic monitoring and 
evaluation of eco-environmental quality in the Three-North region. The Theil-Sen median trend 
method combined with the MK test was used to analyze the variation trend of eco-environmental 
quality in the Three-North region during the study period, and the Hurst exponent and the 
Theil-Sen median trend were superimposed to predict the future trend of eco-environmental 
quality. In addition, the main drivers influencing the spatial distribution of eco-environmental 
quality were identified using a geographical detector to serve as a reference for securing the 
ecological environment in the Three-North region. 


2 Materials and methods 


2.1 Study area 


The study area is the Three-North region (i.e., the three northern regions of China), located 
between 31°31’-53°34'N and 73°26'—-135°06’E (Fig. 1). It includes Northwest China, North 
China and Northeast China, covering 13 provinces, autonomous regions and municipalities, with 
a total area of approximately 5,470,000 km’. The Three-North region encompasses vast 
deserts, sandy lands and the Gobi, most of which are located in arid and semi-arid climate 
zones. The natural environment is fragile in this region. The western part is dominated by 
desert lands, while the eastern part is mainly covered by forestlands and grasslands. Since the 
implementation of the Three-North Shelterbelt Program, a total afforestation and preservation 
area of 301,400 km? has been completed in the project area (Ji et al., 2022). Over the past 40 a, 
the forest coverage has increased from 5.05% to 13.57%, significantly improving the 
ecological conditions (Li and Feng, 2021). 


2.2 Data sources and processing 


The remote sensing data used to calculate the RSEI include the MODO09A1 V6 and MOD11A2 V6 
datasets. MODO9A1 V6 contains the surface reflectance data from the Terra MODIS sensor with 
a spatial resolution of 0.5 km, which is corrected for atmospheric conditions (Xu et al., 2019). 
MODI11A2 V6 contains the surface temperature data with a spatial resolution of 1.0 km. To 
reduce the interference caused by the difference in image acquisition time, we obtained both 
datasets (MOD09A1 V6 and MOD11A2 V6) for the growing period (July-September) from 2000 
to 2020 on the GEE platform. In addition, we declouded the surface reflectance data using the 
quality assessment band to improve the image quality under cloud interference (Zheng and Zhu, 
2017) and synthesized using the median value to stabilize the MODIS images for subsequent 
analysis. 

Changes in ecosystems are influenced by a combination of natural factors and human activities 
(Jiang et al., 2021). Temperature, precipitation and other important meteorological factors have 
multiple impacts on ecosystems (Yuan et al., 2021). Topographic factors, as important 
environmental factors, can indirectly affect the growth and distribution of vegetation by 
influencing climate and soil (Zhang et al., 2021b). Population and socioeconomic development 
are important anthropogenic factors, and land use type and soil health can reflect the state of the 
ecological environment (Yuan et al., 2021). To explore the drivers of the spatial differentiation in 
eco-environmental quality in the study area, we examined ten drivers from two categories of 
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natural and anthropogenic factors. The natural factors included meteorological factors 
(cumulative precipitation, average temperature, average relative humidity and potential 
evapotranspiration) and topographic factors (elevation, slope and aspect). To obtain 
meteorological data for the corresponding time periods, we summarized the monthly average 
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Fig. 1 Location of the study area (Three-North region) in China (a) and topography of the study area (b). The 
Three-North region includes Northwest China, North China and Northeast China, covering 13 provinces, 
autonomous regions and municipalities. SAR, special administrative region. 
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temperature, monthly average relative humidity and monthly potential evapotranspiration data 
from July to September by the mean values, and calculated the monthly cumulative precipitation 
data by the total amount. Anthropogenic factors included socioeconomic factors (nighttime-light 
and population density) and human disturbance factor (land use type). The land use type data 
were classified using a two-level classification system: the first level was divided into six types 
according to land resource and its use attribute; and the second level was divided into 25 types 
according to the natural attribute of land resource (Liu et al., 2010). The six types in the first level 
classification were used for analysis in this study, which included cultivated land, forestland, 
grassland, water bodies, construction land and unused land. More detailed information about 
these data is shown in Table 1. 


Table 1 Detailed description of the data used in the study 
Spatial Time 


Function Data name : : Time span Source 
resolution resolution 
RSEI MODO9A1 V6 0.5 km 8d 2000-2020 GEE 
calculation MODI11A2 V6 1.0km 8d (growing period) (https://developers.google.com) 
Cumulative precipitation 1.0 km Monthly National Earth System Science 
Average temperature 1.0km Monthly 2000, 2005, 2010, Data Center and the National 
: as 2015 and 2020 Science & Technology 
Average relative humidity 1.0km Monthly (growing period) Infra-structure of China 
Potential evapotranspiration 1.0 km Monthly (http://www.geodata.cn) 
Eevatión 1-0 km 2000 A Big Earth Data Platform for 
Slope 1.0 km 2000 Three Poles 
Factor (http://poles.tpdc.ac.cn/) 
response Aspect 1.0 km 2000 
analysis A Big Earth Data Platform for 
Nighttime-light 1.0 km Annual e pee i 0, Three Poles 
(http://poles.tpdc.ac.cn/) 
. A 2000, 2005, 2010, WorldPop 
Population density Loki G 2015 and 2020 (https://www.worldpop.org) 
Resource and Environment 
Land use type 1.0km Quinquennial 2000, 2005, 2010, Science and Data Center of the 


2015 and 2020 Chinese Academy of Sciences 
(https://www.resde.cn) 


Note: RSEI, remote sensing ecological index; GEE, Google Earth Engine. Growing period refers to months from July to 
September. 


2.3 Methods 


2.3.1 RSEI construction 


The RSEI can assess eco-environmental quality by greenness, wetness, dryness and heat (Hu and 
Xu, 2018). In this study, the NDVI was used as the greenness indicator; the humidity component 
(Wet) obtained through the tasseled cap transformation was used as the wetness indicator (Xu et al., 
2021); the normalized difference built-up and soil index (NDBSJ), i.e., the mean of the soil index 
(SI) and index-based built-up index (IBI), was used as the dryness indicator; and the LST was used 
as the heat indicator. 

On the GEE platform, the NDVI, Wet and NDBSI parameters were extracted using the processed 
surface reflectance data of MODIS images and resampled to a spatial resolution of 1.0 km to 
maintain consistency between different data sources. The calculation formulas of each index are as 
follows: 


NDVI =P., (1) 
P2+ Pi 


Wet = 0.3279p, +0.3406p, +0.1509/, + 0.1973, — 0.7112, —0.4572p,, (2) 
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SI ; (3) 
(ps +1) +(P2 +p) 
206 Pr P4 l 
pr- e tA Pate , (4) 
2 P6 j l P2 P4 l 
Pe Pr \P2tPi Pat Po 
NDBSI = SE (5) 


where NDVI is the normalized difference vegetation index; Wet is the humidity component; pi, 
p2, P3, P4, Po and p7 are the reflectance values of the red, near-infrared 1, blue, green, short 
wavelength infrared 1 and short wavelength infrared 2 bands of MODO9A1, respectively; SI is the 
soil index; IBI is the index-based built-up index; and NDBSI is the normalized difference built-up 
and soil index. 

The daytime surface temperature data of MOD11A2 V6 were synthesized to obtain the LST. 
All indicators were standardized to eliminate the influence of the dimension between different 
indicators. In addition, water bodies were masked using the JRC Yearly Water Classification 
History dataset (Pekel et al., 2016) provided by the GEE platform (Source: EC JRC/Google), and 
the union of different years of water bodies within the study area was masked in the 
spatiotemporal analysis. Principal component analysis (PCA) was used to integrate the four 
indicators (Hong et al., 2020), and the first principal component was defined as the RSEIo. The 
RSEIo was then standardized to obtain the RSEI to facilitate the measurement and comparison of 
the four indicators, as follows: 


RSEI, = f (NDVI Wet, NDBSLLST), (6) 
RSET = RSElo -RSElomin a 
RSET max = RSE] omin 


where RSEIo is the first principal component of the four indicators; LST is the land surface 
temperature (K); RSEI is the remote sensing ecological index; RSElomin is the minimum value of 
the RSEIo; and RSEIomax is the maximum value of the RSEIo. 

The RSEI value is between 0.0 and 1.0. The closer the value is to 1.0, the better the quality of 
the ecological environment; and the closer the value is to 0.0, the worse the quality of the 
ecological environment. The RSEI values in the study were divided into five intervals (0.0—0.2, 
0.2—-0.4, 0.40.6, 0.6-0.8 and 0.8—1.0), corresponding to the five grades of eco-environmental 
quality: poor, fair, moderate, good and excellent, respectively. 

2.3.2 Trend analysis 


The Theil-Sen median trend method combined with the MK test was used to analyze the variation 
trend of eco-environmental quality in the Three-North region during the study period, and the 
Hurst exponent and the Theil-Sen median trend were superimposed to predict the future trend of 
eco-environmental quality. 
(1) Trend analysis based on the Theil-Sen median trend method and MK test 

As a robust method of trend calculation, the Theil-Sen median trend method is often used in 
conjunction with the MK test to detect the variation trend of long-term sequence data (Li et al., 
2021a; Jiang et al., 2022). The calculation formula of the Theil-Sen median (£) is as follows: 


RSEI , — RSEI; 


B= median[ ) (2000 <i < j < 2020), (8) 


where i and j represent the RSEI time series (i.e., 21 a from 2000 to 2020); and RSEI; and RSEĻ; 
represent the RSEI values of years i and j, respectively. 
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Referring to the classification standard of Jiang et al. (2015), we divided the variation trend of 
eco-environmental quality into three types based on the Theil-Sen median (£): f>0.0005 
indicating an improved trend of eco-environmental quality; —0.0005<$<0.0005 indicating a stable 
trend of eco-environmental quality; and < —0.0005 indicating a degraded trend of eco- 
environmental quality. The MK method used the Z statistic for the significance test. Based on a 
significance level of a, when |Z|>Zı-a2, the trend of the time series is considered significant; 
otherwise, it is considered nonsignificant. Referring to Li et al. (2021b), the value of a was taken 
as 0.05 in this study, which indicated that the trend is judged to be significant at a confidence 
level of 95%. When a equals to 0.05, the Z-score corresponding to the cumulative probability of 
1—a/2 is queried through the standard normal distribution table as 1.96 (Kutner et al., 2004). That 
is, the value of Zı-a⁄2 was 1.96 in this study. Combining the Theil-Sen median and MK test results, 
we divided the variation trend of eco-environmental quality into five categories. The specific 
classification standard is shown in Table 2. 


Table 2 Classification standard for the trend analysis of eco-environmental quality based on the Theil-Sen 
median (£) and Mann-Kendall (MK) test results 


Classification standard Variation trend of eco-environmental quality 
£< -0.0005 and |Z|>1.96 Significant degradation 
2< 0.0005 and |Z|<1.96 Slight degradation 
—0.0005<£<0.0005 Basically stable 
f>0.0005 and |Z|<1.96 Slight improvement 
f>0.0005 and |Z|>1.96 Significant improvement 


Note: The MK method used the Z statistic for the significance test. 


(2) Hurst exponent 

The Hurst exponent reflects the autocorrelation of time series and can effectively predict future 
development trend (Li et al., 2021b; Rivas-Tabares et al., 2021). The calculation procedure is 
shown as follows. 

Given a time sequence RSEI(f) (where 1, 2, ..., n), the following equations can be obtained 
for any positive integer (t21, where T is the integer). 


Mean sequence: (RSEI)_ = iy RSEI(t) (t =1, 2, ..., n). (9) 
T = 
Accumulated deviation: X(t,7)= > (RSEI(¢)—(RSEI),) (<t 7). (10) 
t= 
Range: R(T)= max X(t,t)- min X(¢,r) (T=1, 2, ..., n). (11) 
Isisr I<t<t 
1< Vy 
Standard deviation: S(r)= -5 (RSEI(¢) — (RSEI), y (t=1, 2, ..., n). (12) 


t=l 

If there exists R(1)/S(t)°= 1”, then the Hurst phenomenon exists in time series (RSEI(¢)), where 
H is the Hurst exponent. The Hurst exponent value can be fitted by the least square method in 
double logarithmic coordinate system (In(z), In(R(z)/S(z))), and its value ranges from 0.00 to 1.00. 

When 0.00<H<0.50 (where H is the value of Hurst exponent), the long-term correlation of 
eco-environmental quality in time series is characterized by anti-sustainability; the smaller the 
value, the stronger the anti-sustainability. When H=0.50, there is no long-term correlation of 
eco-environmental quality in time series. When 0.50<H<1.00, the long-term correlation in time 
series is characterized by persistence; the larger the value, the stronger the persistence. The Hurst 
exponent and the Theil-Sen median trend were superimposed to analyze the future trend of 
eco-environmental quality. We classified the future trend of eco-environmental quality following 
the classification standard of Zhang et al. (2021b), as shown in Table 3. 
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Table 3 Classification standard for the future trend analysis of eco-environmental quality based on the Hurst 
exponent and the Theil-Sen median trend 


Classification standard Future trend of eco-environmental quality 
£>0.0005 and 0.65<H<1.00 Strong persistent improvement 
f>0.0005 and 0.50<H<0.65 Weak persistent improvement 

£< —0.0005 and 0.00<H<0.35 Strong anti-sustained degradation 
f<-0.0005 and 0.35<H<0.50 Weak anti-sustained degradation 
—0.0005<£<0.0005 Essentially constant 
$< 0.0005 and 0.65<H<1.00 Strong persistent degradation 
$< -0.0005 and 0.50<H<0.65 Weak persistent degradation 
£>0.0005 and 0.00<H<0.35 Strong anti-sustained improvement 
f>0.0005 and 0.35<H<0.50 Weak anti-sustained improvement 


Note: H is the value of Hurst exponent. 


2.3.3 Geographical detector 


The geographical detector, which can measure the spatial differentiation and reveal the driving 
force behind the differentiation mechanism, has been widely used to analyze the influencing 
mechanisms of various factors (Wang and Hu, 2012; Su et al., 2020). In this study, the RSEI was 
selected as the dependent variable, and cumulative precipitation, average temperature, average 
relative humidity, potential evapotranspiration, elevation, slope, aspect, nighttime-light, 
population density and land use type were selected as the independent variables. First, a 5.0 
kmx5.0 km fishing network was used to uniformly generate 207,034 points in the Three-North 
region, through which the dependent variable was matched with the independent variables. Then, 
each independent variable in 2000, 2005, 2010, 2015 and 2020 was converted into a type quantity, 
and each hierarchical interval was determined by the optimal q statistic (where q is the 
explanatory power). Among them, as a type quantity, the six types in the first level classification 
were used for analysis in this study. Finally, the effects of each factor on the spatial distribution of 
the RSEI in 2000, 2005, 2010, 2015 and 2020 were investigated by factor detector, interaction 
detector and ecological detector of the geographical detector. 

The factor detector evaluated the effect of different environmental factors on the spatial 
differentiation of eco-environmental quality by calculating the q value. The range of the q value is 
from 0.000 to 1.000; the higher the g value, the greater the influence. 


L 
q _ 1 z ae N,On 
No? 


; (13) 


where q represents the explanatory power; L is the number of layers of the factor; h is the 
classification or partition of factor Y or factor X; N, and N are the unit numbers of class A and the 
entire region, respectively; and a,” and o? are the variances of Y for class A and the entire region, 
respectively. 

The interaction detector identified interaction between two factors to evaluate whether these 
two factors affect eco-environmental quality jointly or independently (Table 4). 


Table 4 Classification standard of the interaction type of factors 


Classification standard Interaction type 
q(X1NX2)<Min[g(X1), g(X2)] Nonlinear-weaken 
Min[q(X1), g(X2)]<q(X1N.X2)<Max[q(X1), g(X2)] Uni-weaken 
q(X1NX2)>Max[q(X1), q(X2)] Bi-enhance 
q(X1NX2)=q(X1)+q(X2) Independent 
Q(X1N.X2)>q(X1)+q(X2) Nonlinear-enhance 


Note: g(X1) and q(X2) represent the explanatory power (q value) of factors X1 and X2 for eco-environmental quality, 
respectively; g(X1MX2) represents the explanatory power (q value) for eco-environmental quality when factors X1 and X2 
interact. 
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The ecological detector was used to determine whether the two factors had significant 
differences in the spatial distribution of eco-environmental quality. This was tested using the F 
statistic as follows: 


Ll 
F= Nyi (Nx2 -DX a Nor 
E L2 á 
N y2(N xı -DX Nor 


where Nx, and Ny represent the sample size of the two factors, respectively; and L1 and L2 
represent the number of the layers of the two factors, respectively. 


(14) 


3 Results 


3.1 Spatial distribution of eco-environmental quality 


Figure 2 shows the spatial distribution of eco-environmental quality in the study area for some 
typical years from 2000 to 2020. The RSEI values in Northwest China were generally between 
0.2 and 0.4; it had been at a low level for a long time. The eco-environmental quality in the 
northern slope valley of the Tianshan Mountains and the Qaidam Basin was relatively good, and 
the RSEI values were mostly between 0.4 and 0.8. Most of these areas had a certain degree of soil 
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Fig. 2 Spatial distribution of eco-environmental quality (as indicated by the RSED) in the Three-North region in 
2000 (a), 2002 (b), 2004 (c), 2006 (d), 2008 (e), 2010 (£), 2012 (g), 2014 (h), 2016 (i), 2018 (j) and 2020 (k) 
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moisture, and the surface temperature was lower than that in the desert area. The RSEI values in 
North China gradually increased from north to south. The RSEI values in the Alxa Plateau, the 
northern Hetao Plain and the Mu Us Sandy Land were mainly between 0.2 and 0.4, and those in 
the Hulun Buir Plateau and Otindag Sandy Land were mainly in the range of 0.4—0.6. The 
eco-environmental quality of the Haihe River Basin and its surrounding areas in the southern part 
of North China was mainly classified as good (the RSEI values ranging from 0.6 to 0.8). The 
eco-environmental quality of Northeast China was mostly classified as good (the RSEI values of 
0.6-0.8) or excellent (the RSEI values of 0.8-1.0). The RSEI values of the Yilehuli Mountain, 
Lesser Khingan Mountains and Changbai Mountains were relatively high, generally higher than 
0.8. 

The RSEI values in the Three-North region from 2000 to 2020 were classified and counted 
(Fig. 3). There were obvious differences in the percentage of the area with various RSEI values in 
different regions. Approximately 39.00%-45.50% of the area in Northwest China had RSEI 
values between 0.2 and 0.4, while approximately 21.80%-—33.20% of the area had RSEI values 
greater than 0.6. Compared with Northwest China, the percentage of the area in North China with 
RSEI values below 0.4 was smaller while that with the RSEI values between 0.6 and 0.8 was 
greater. The RSEI values in Northeast China were all greater than 0.2. In Northeast China, the 
area with the RSEI values ranging from 0.2 to 0.4 accounted for approximately 0.20%—0.50% of 
the total area in this region, in which eco-environmental quality was considered as fair, while 
eco-environmental quality in the remainder region was classified as moderate, good or excellent. 
From 2000 to 2020, the percentage of the area classified at each RSEI level in the study area 
changed to different degrees, but the RSEI exhibited a fluctuating upward trend overall and its 
annual mean growth rate was 0.0016/a. The mean RSEI values were between 0.5 and 0.6 from 
2000 to 2020 in the Three-North region, with the minimum value occurring in 2000 and the 
maximum in 2012. 


Eco-environmental quality 


= Poor (0.0<RSEI<0.2) = Fair (0.2<RSEI<0.4) = Moderate (0.4<RSEI<0.6) 
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Fig. 3 Percentage of the area with different eco-environmental quality degrees (as indicated by the RSEI values) 
and the mean RSEI in Northwest China (a), North China (b), Northeast China (c) and the Three-North region (d) 
from 2000 to 2020 
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3.2 Dynamic tendency of eco-environmental quality 


The statistical results of the variation trend of eco-environmental quality from 2000 to 2020 are 
shown in Table 5. During the study period, approximately 56.92% of the area in Northwest China 
showed an improving trend in eco-environmental quality. According to the RSEI, Northwest 
China also contained the largest area with slight improvement in eco-environmental quality, at 
1,051,168 km? (accounting for 34.94% of Northwest China), followed by the area with basically 
stable eco-environmental quality, accounting for 34.13% of Northwest China. More than 75.00% 
of the areas in North China and Northeast China respectively showed an improving trend in 
eco-environmental quality. The areas with slightly and significantly improved eco-environmental 
quality in North China were 677,542 km? and 491,153 km?, accounting for 44.70% and 32.41% of 
the total area in this region, respectively. In Northeast China, the areas with slight and significant 
improvement in eco-environmental quality were 256,214 km? and 335,402 km?, respectively, at 
percentages of 32.55% and 42.61%, respectively. The significantly degraded areas in Northwest 
China, North China and Northeast China were 41,949, 21,793 and 11,238 km’, respectively, 
accounting for 1.41% of the total study area (the Three-North region). In terms of the overall 
values, 7.79% of the total area exhibited a degradation trend in eco-environmental quality in the 
Three-North region, and the main trend of eco-environmental quality from 2000 to 2020 was 
basically stable or improved (accounting for 26.82% and 65.39%, respectively). 


Table 5 Statistical results of the variation trend in eco-environmental quality from 2000 to 2020 


Northwest China North China Northeast China Three-North region 
Variation trend Area Percentage Area Percentage Area Percentage Area Percentage 
(km’) (%) (km?) (%) (km?) %) (km’) (*) 
Significant degradation 41,949 1.39 21,793 1.44 11,238 1.43 75,190 1.41 
Slight degradation 227,387 7.56 78,656 5.19 32,731 4.16 339,363 6.38 
Basically stable 1,026,679 34.13 246,460 16.26 151,471 19.25 1,426,634 26.82 


Slight improvement 1,051,168 34.94 677,542 44.70 256,214 32.55 1,987,808 37.38 
Significant improvement 661,078 21.98 491,153 32.41 335,402 42.61 1,489,976 28.01 


According to the spatial distribution of the variation trend of eco-environmental quality (Fig. 
4), the areas exhibiting a basically stable trend of eco-environmental quality in Northwest China 
were distributed in the Tarim Basin, Junggar Basin, Kumtag Desert and Qaidam Basin. Similarly, 
eco-environmental quality in the Badain Jaran Desert, Alxa Plateau and the surrounding areas of 
Central Gobi also remained basically stable. These areas are mainly desert and Gobi regions with 
low soil moisture content, high surface temperature, and poor eco-environmental quality for a 
long time. With high vegetation coverage and relatively stable ecosystems, eco-environmental 
quality of the Lesser Khingan Mountains and Changbai Mountains in Northeast China remained 
basically stable at a good level. The areas north of the Yarkant River Basin and Tarim River 
Basin, the Loess Plateau and the Northeast China Plain (Songnen Plain, Liaohe River Plain and 
Sanjiang Plain) mainly showed a significant improvement in eco-environmental quality. However, 
there were some areas showing significant degradation trend of eco-environmental quality in 
Northwest China, North China and Northeast China. In Northwest China, these areas were mainly 
distributed in the Halktawu Mountain and Boluokenu Mountain (Fig. 4a); in North China, they 
were mainly distributed in areas with rapid economic development, such as Beijing-Tianjin-Hebei 
region (Fig. 4b). In addition, the northernmost end of the border between Inner Mongolia 
Autonomous Region and Heilongjiang Province (Fig. 4c), that is, the eastern side of the Ergun 
River and the western side of the Mohe City, also showed significant degradation in 
eco-environmental quality. 

The Hurst exponent was superimposed and analyzed with the results of Theil-Sen median trend 
analysis to obtain the persistence of the future trend in eco-environmental quality in the 
Three-North region (Fig. 5). The Hurst exponent of the RSEI in the Three-North region ranged 
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Fig. 4 Spatial distribution of the variation trend of eco-environmental quality in the Three-North region from 
2000 to 2020. The pie chart shows the percentage of the area occupied by different variation trends of 
eco-environmental quality. (a), (b) and (c) show the partial enlargement of the spatial distribution of the variation 
trend of eco-environmental quality. 


from 0.20 to 0.76, with an average of 0.53, and the percentages of the area with persistence and 
anti-sustainability in eco-environmental quality were 63.55% and 36.44%, respectively (Fig. 5a). 
The overall future trend of eco-environmental quality in the Three-North region was weak, 
mainly dominated by weak persistence, followed by weak anti-sustainability, accounting for 
52.47% and 35.24% of the total study area, respectively, indicating that the future trend of 
eco-environmental quality in most areas may be reversed. The percentage of the area with strong 
persistence in eco-environmental quality was 11.08%, mainly located in the eastern Loess 
Plateau, the Songnen Plain and the Liaohe River Plain. As shown in Figure 5b, 0.81% of the 
Three-North region showed strong persistent degradation in eco-environmental quality, 0.31% 
exhibited strong anti-sustained improvement, 10.09% showed strong persistent improvement and 
0.13% showed strong anti-sustained degradation. The overall future changes of 
eco-environmental quality in the Three-North region were mainly based on improvement in the 
RSEI, but they were not strong in sustainability. The percentages of the area with weak 
anti-sustained improvement and weak persistent degradation in eco-environmental quality (the 
RSEI might decrease) in the Three-North region were 16.23% and 3.78%, respectively. The areas 
with weak persistent improvement and weak anti-sustained degradation in eco-environmental 
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Fig. 5 Spatial distribution of the long-term correlation of eco-environmental quality in time series (indicated by 


the Hurst exponent (H); a) and the future trend of eco-environmental quality (b) in the Three-North region 
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quality (the RSEI might increase) accounted for 39.06% and 2.75% of the Three-North region, 
respectively. In the future, the areas where the RSEI would remain essentially constant are mainly 
distributed in the Tarim Basin, the Junggar Basin, the Qaidam Basin, the Badain Jaran Desert, the 
Alxa Plateau and the surroundings of the Central Gobi; these expectations are highly consistent 
with the variation trend over the past 21 years (2000-2020). The eco-environmental quality of the 
Halktawu Mountain, Boluokenu Mountain and Beijing-Tianjin-Hebei region may also continue to 
degrade to varying degrees in the future. 


3.3 Analysis between influencing factors and eco-environment quality 


Ten influencing factors, including cumulative precipitation, average temperature, average relative 
humidity, potential evapotranspiration, elevation, slope, aspect, nighttime-light, population 
density and land use type, were selected and their impacts on the spatial distribution of 
eco-environmental quality in the Three-North region were explored by the geographical detector. 
The explanatory power q values of the influencing factors were obtained, as shown in Figure 6. 
Based on these results, the driving effect of each influencing factor on the spatial distribution 
pattern of eco-environmental quality was different. From 2000 to 2020, the q value of average 
relative humidity was the largest among the ten influencing factors, followed by cumulative 
precipitation, land use type and potential evapotranspiration. The q value of average relative 
humidity was greater than 0.870, indicating that average relative humidity had a great impact on 
the formation and change in the spatial distribution of eco-environmental quality. In addition, 
cumulative precipitation and land use type were the subdominant factors affecting the spatial 
distribution of eco-environmental quality in the Three-North region, with q values of 
approximately 0.730 and 0.560, respectively. The g values of potential evapotranspiration 
fluctuated around 0.410 and decreased to a greater extent in 2015. The other influencing factors 
had relatively small effects on the spatial distribution of eco-environmental quality, with q values 
less than 0.300. Among them, aspect had the weakest effect on the spatial distribution pattern of 
eco-environmental quality, with q value of only approximately 0.010. In general, the explanatory 
power of meteorological factors on the spatial distribution of eco-environmental quality was 
stronger than that of topographic and socioeconomic factors. The effects of anthropogenic factors, 
such as population density and land use type, on eco-environmental quality gradually increased 
over time. 

A two-factor interaction would enhance the explanatory power of influencing factors on the 
spatial distribution of eco-environmental quality (Fig. 7). From 2000 to 2020, the interaction 
types between the influencing factors were mainly double synergy (bi-enhance) and nonlinear 
synergy (nonlinear enhance), and there were no factors to act independently. We arranged the 
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Fig. 6 Explanatory power q values of cumulative precipitation (X1), average temperature (X2), average relative 
humidity (X3), potential evapotranspiration (X4), elevation (X5), slope (Xo), aspect (X7), nighttime-light (Xs), 
population density (Xo), and land use type (X10) in 2000, 2005, 2010, 2015 and 2020 
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Fig. 7 Interaction detection and ecological detection results of cumulative precipitation (X1), average 
temperature (X2), average relative humidity (X3), potential evapotranspiration (X4), elevation (X5), slope (X6), 
aspect (X7), nighttime-light (Xs), population density (Xo) and land use type (X10) in 2000 (a), 2005 (b), 2010 (c), 
2015 (d) and 2020 (e). The f means that the interaction between two influencing factors is nonlinear enhance 
effect and no f means that the interaction between two influencing factors is bi-enhance effect. 


explanatory force of the factor interaction in five representative years according to the g values. 
The most dominant interaction factors were the combination of average temperature and average 
relative humidity in 2000, 2005, 2010, 2015 and 2020, with q values of 0.920, 0.922, 0.918, 0.923 
and 0.929, respectively. In addition, the q values of average relative humidityMpotential 
evapotranspiration, average relative humidityNelevation, and average relative humidityNland use 
type were all greater than 0.900. Average relative humidity was the dominant meteorological 
factor on the spatial distribution of eco-environmental quality and made the greatest contribution 
to the change in the distribution pattern, including interactions with other influencing factors. The 
q values of elevation and population density were low, but the explanatory power between them 
was stronger than that of the other single factors when interacting with meteorological factors. 
The interaction between slope and other influencing factors was mostly nonlinearly enhanced. In 
addition, the explanatory power of topographic and socioeconomic factors was weak when acting 
as single factors, and the power remained low after a two-factor interaction. Topographic and 
socioeconomic factors presented bi-enhance effect and nonlinear enhance effect when interacting 
with meteorological factors, indicating that the influence of topographic and socioeconomic 
factors on the spatial distribution of eco-environmental quality was neither unilateral nor simply 
superposed with meteorological factors. In addition, the ecological detection of each factor was 
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performed to detect whether there were significant differences in the impact of the two factors on 
the spatial distribution of eco-environmental quality (Fig. 7). According to the detection results, 
except for average temperature and elevation that were not significantly different in 2015, all the 
other factors differed significantly and had different influence mechanisms on the spatial 
distribution of eco-environmental quality. 


4 Discussion 


4.1 Variation trend of eco-environmental quality 


According to the variation trend analysis results (Fig. 4), areas with significant degradation in 
eco-environmental quality were found in Northwest China, North China and Northeast China 
from 2000 to 2020. In Northwest China, these areas were mainly concentrated on both sides of 
the geometric center of the Tianshan Mountains. Zhou et al. (2020) showed that high altitude, low 
precipitation, dry climate and fragile vegetation are the main causes for the degradation of 
eco-environmental quality in this area. Meanwhile, there were obvious regional differences in 
eco-environmental quality in the Beijing-Tianjin-Hebei region of North China, generally 
exhibiting an improvement pattern to the northwest while a degradation pattern to the southeast. 
Due to the strengthening of artificial ecological protection policies, eco-environmental quality in 
Zhangjiakou City, Chengde City and other places has been improved (Deng et al., 2018; Ji et al., 
2020). Although ecological restoration projects are conducive to the benign development of 
eco-environmental quality (Xu et al., 2020), areas with dense populations, convenient 
transportation and strong human disturbance are not suitable for vegetation restoration (Zhang et 
al., 2017), resulting in a decline in eco-environmental quality in the southeastern Hebei Province 
(such as Handan and Hengshui cities) and Tianjin. Located in the northern end of China, a small 
area on the east side of the Ergun River experienced a decline in eco-environmental quality due to 
extreme climatic conditions. Since the implementation of the Three-North Shelterbelt Program, 
the overall eco-environmental quality of the Three-North region was improved. However, a 
significant degradation trend in eco-environmental quality in some regions existed, and 
eco-environmental quality might continue to degrade in the future (Fig. 5b). The change in 
eco-environmental quality was affected by many factors, such as climate change, land use and 
population density (Yuan et al., 2021), but the reasons for the degradation in eco-environmental 
quality varied from region to region. In the future, ecological protection policies should be 
formulated according to the main factors affecting the degradation of eco-environmental quality 
in each area. For areas with harsh climatic conditions and fragile vegetation ecology, such as the 
edge of the oases in the northern and southern Tianshan Mountains, attention should be paid to 
ecological water use and strict protection of the natural surroundings (Zheng and Zhu, 2017). For 
areas with rapid economic development, such as the Beijing-Tianjin-Hebei region, some measures 
should be considered to restore key vegetation in regions with intensive human disturbance (Yuan 
et al., 2021). For areas with better natural conditions and abundant resources, such as the Greater 
Khingan Mountains and Changbai Mountains, priority should be given to protection and 
categorical management to promote ecological restoration and functional reconstruction (Wang et 
al., 2020). 


4.2 Factors influencing the spatial distribution of eco-environmental quality 


The interaction between natural factors is complex, and it is difficult to ensure that each factor is 
independent of each other. Therefore, analytical methods will be considered relatively more in the 
driver analysis to reveal the independent effects of factors as much as possible. In contrast to 
correlation and partial correlation analyses, the geographical detector has no assumption of 
linearity and is essentially free of multicollinearity, with each factor having an independent effect 
on the result. The geographical detector can diagnose the independent driving force of 
geographical variables, but its application in eco-environmental quality analysis is still relatively 
rare. In this study, ten influencing variables from two categories of natural and anthropogenic 
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factors were analyzed to determine the drivers of the spatial differentiation of eco-environmental 
quality by the geographical detector. To make the result of the geographical detector as accurate 
as possible, we uniformly generated a total of 207,034 points in the Three-North region to match 
the RSEI with each influencing factor. The hierarchical interval of continuous variables has direct 
or indirect effects on the q value (Wang et al., 2010; Su et al., 2020). Therefore, in the 
discretization of the continuous variables, the q statistic was used to prefer parameters, such as 
discrete methods and the number of levels, and determine the hierarchical interval of each 
influencing factor. The larger the q value, the better the partition effect (Su et al., 2020). 

Based on the factor detection results (Fig. 6), average relative humidity and accumulated 
precipitation had strong explanatory power for the spatial distribution of eco-environmental 
quality, while average temperature had weak explanatory power. Most of the Three-North region 
is located in arid and semi-arid climate zones. Temperature mainly regulated the annual growth 
law of vegetation. Compared with precipitation, the correlation between temperature and 
interannual variation in vegetation was generally low, and vegetation was more sensitive to the 
influence of precipitation and humidity (Xie et al., 2016). In addition, the influence of 
anthropogenic factors such as population density and land use type on the spatial distribution of 
eco-environmental quality gradually increased over time, which is consistent with the findings of 
Chun et al. (2018) and Ji et al. (2020), supporting the reasonableness of the factor detection 
results. The results of interaction detection (Fig. 7) showed that the interaction between each two 
drivers (factors) had a greater impact on the spatial distribution of eco-environmental quality than 
a single driver (factor). Topographic factors (e.g., elevation and slope) had a moderate effect on 
the spatial distribution of eco-environmental quality under the action of a single factor, but their 
degree of influence was enhanced when interacting with meteorological factors. This suggested 
that multiple drivers (factors) comprehensively affect the quality of the regional eco-environment 
(Dai et al., 2020). However, when exploring the explanatory power of each influencing factor on 
the spatial distribution of eco-environmental quality, we only discussed the global factor 
responses of five key years (2000, 2005, 2010, 2015 and 2020) in the entire Three-North region 
without further dividing the study area or analyzing the local driving forces of different regions 
on an annual scale. In follow-up studies, the spatial and temporal resolution of the analysis should 
be improved, and the specific response of each influencing factor in different regions should be 
discussed in depth. 


4.3 Limitations and prospects 


GEE is capable of rapidly batch-processing large numbers of MODIS or Landsat images, which is 
an improvement on the tedious calculation process of the RSEI in large-scale applications (Huang 
et al., 2021; Xiong et al., 2021). In this study, the RSEI was used to evaluate eco-environmental 
quality in the Three-North region. The four indicators (NDVI, Wet, NDBSI and LST) used to 
construct the RSEI had certain loadings on the first principal component. Unlike the positive and 
negative uncertainties of the loadings of the four indicators on the second, third and fourth 
principal components, the eigenvectors of the greenness and humidity indicators on the first 
principal component were always positive, while the eigenvectors of dryness and heat were 
always negative (data not shown). In addition, the eigenvalue of the first principal component had 
a contribution rate of more than 70.00% in all years, and it integrated most of the characteristics 
of the four indicators, which could reflect the eco-environment condition more reasonably. The 
spatial distribution results (Fig. 3) showed that the RSEI values in the Three-North region were 
generally moderate, and the percentage of the area with particularly good or poor 
eco-environmental quality was relatively small. The RSEI values were greater in Northeast China 
and lower in Northwest China, which is consistent with the results obtained by Liao et al. (2020) 
and Ji et al. (2022). The RSEI can effectively distinguish the ecological status of different 
landscapes (Wang et al., 2010), showing its effectiveness for evaluating the changes in 
eco-environmental quality. However, it still has some limitations. Due to the large scope of the 
study area, it was difficult to unify the dates of the selected MODIS images when constructing the 
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RSEI. Although the median synthesis method was designed to ensure the intermediate level of 
temporal remote sensing images as much as possible, it might have some impact on the 
comparability of the RSEI evaluation (Zheng et al., 2020). There has yet to be a complete solution 
for the instability of the RSEI time series, which greatly limits its utility in large-scale 
applications (Zheng et al., 2022). In addition, the evaluation method of eco-environmental quality 
based on the RSEI mainly focused on the four indicators of greenness, wetness, dryness and heat. 
However, the ecosystems are complex and diverse, involving many aspects related not only to the 
natural environment but also to the socioeconomic development and government management 
(Xiong et al., 2021; Yuan et al., 2021). How to effectively screen other key indicators and 
combine them with remote sensing technology is a direction for future research. 


5 Conclusions 


In this study, the Three-North region was selected as the research area, and MODIS images 
obtained from the GEE platform were used to calculate the RSEI. The Theil-Sen median trend 
method combined with the MK test was used to analyze the variation trend of eco-environmental 
quality in the Three-North region from 2000 to 2020, and the Hurst exponent and the Theil-Sen 
median trend were superimposed to predict the future trend of eco-environmental quality. In 
addition, the geographical detector was applied to quantitatively analyze the influence of different 
factors on the spatial distribution of eco-environmental quality. 

The spatial distribution of the eco-environmental quality in the Three-North region exhibited 
obvious regional characteristics. Eco-environmental quality had been at a low level for a long 
time in Northwest China, gradually increased from north to south in North China, and was 
relatively good in Northeast China. From 2000 to 2020, the average RSEI values in the 
Three-North region increased at an average growth rate of 0.0016/a. Eco-environmental quality in 
the Three-North region showed the spatial distribution characteristics of overall improvement and 
local degradation. The persistence of future change in eco-environmental quality was stronger 
than the anti-sustainability. Eco-environmental quality may be mainly improved in the future, but 
the overall trend is weak, and there is a possibility of reversal in some local regions. Thus, it is 
necessary to continue to strengthen the protection of the ecological environment in the 
Three-North region. Average relative humidity, accumulated precipitation and land use type were 
the dominant factors determining the spatial distribution of eco-environmental quality in the 
Three-North region. The explanatory power of meteorological factors on the spatial distribution 
of eco-environmental quality was stronger than that of topographic factors. With passage of time, 
the influence of anthropogenic factors such as population density and land use type on 
eco-environmental quality gradually increased. The two-factor interaction was dominated by the 
bi-enhance effect, and the interaction strengthened the explanatory power of the spatial 
distribution of eco-environmental quality, among which the interaction of average 
temperatureNaverage relative humidity was the highest. 

In the Three-North region, the ecological restoration project is conducive to the benign 
development of eco-environmental quality; however, there are still some areas with degraded 
eco-environmental quality. Generally, eco-environmental quality is influenced by a combination 
of driving factors, and in the future, ecological protection policies should be formulated according 
to the main factors affecting the degradation of eco-environmental quality in various areas. This 
study reasonably reflects the eco-environmental quality of the Three-North region by constructing 
the RSEI, but the stability of the RSEI time series needs further improvement for large-scale 
applications. 
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